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ABSTRACT 

We present a new semi-empirical relation for the mass loss of cool stellar 
winds, which so far has frequently been described by "Reimers' law" . Originally, 
this relation was based solely on dimensional scaling arguments without any 
physical interpretation. In our approach, the wind is assumed to result from the 
spill-over of the extended chromosphere, possibly associated with the action of 
waves, especially Alfven waves, which are used as guidance in the derivation of the 
new formula. We obtain a relation akin to the original Reimers law, but which 
includes two new factors. They reflect how the chromospheric height depends 
on gravity and how the mechanical energy flux depends, mainly, on effective 
temperature. The new relation is tested and sensitively calibrated by modelling 
the blue end of the Horizontal Branch of globular clusters. The most significant 
difference from mass loss rates predicted by the Reimers relation is an increase 
by up to a factor of 3 for luminous late-type (super-)giants, in good agreement 
with observations. 



Subject headings: stars: chromospheres — stars: late-type — stars: mass loss - 
turbulence — waves 
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1. Introduction 

Empirical mass loss formulas are pivotal for the construction of empirical and semiem- 
pirical stellar atmosphere and wind models, stellar evolution computations and studies of 
the interstellar medium, among other topics. Historically, the mass loss rate M of late- 
type giants and supergiants has been described by "Reimers' law", given as M = 77 • ^j^- 
(Reimers 1975, 1977), with L*, R*, M* as stellar luminosity, radius, and mass, respectively, 
given in solar units, and r\ is a fitting parameter. Other empirical mass loss formulas have 
been presented by Lamers (1981), de Jager, Nieuwenhuijzen, & van der Hucht (1988), and 
Nieuwenhuijzen & de Jager (1990), but they do not distinguish between the strong, and now 
well-described dust-driven winds (e.g., Wachter et al. 2002), and the physically very different 
case of nondust- driven winds. 

Despite its wide-ranging success, the mass loss formula by Reimers suffers from two 
important deficiencies. First, it is solely based on dimensional scaling arguments without 
any physical interpretation. In particular, the appearance of the stellar luminosity in the 
formula is awkward noting that for cool star winds, with the exception of molecule-driven 
and dust-driven winds, the luminosity of the star is not expected to be relevant (e.g., Holzer 
& MacGregor 1985). In fact, the Reimers law seems to suggest that a certain fraction of 
the stellar luminosity L* is utilized to lift the wind material from the photosphere. The 
second deficiency consists in the necessity of adjusting the fitting parameter 77, as different 77 
values are required ad hoc to match observed mass-loss rates from different types of giants 
and supergiants. The same is true, if reasonable mass loss yields and final masses are to 
be achieved through stellar evolution models with prescribed mass loss. For more evolved 
AGB giants, the Reimers relation is better replaced by, e.g., de Jager et al. (1988), which 
suggests up to three times as much mass-loss for the tip- AGB (Schroder & Sedlmayr 2001). 
More recently, the Reimers relation also fails to describe revised mass loss rates from K and 
M giant stars, based on updated Ca II ionization balances which consider photoionization 
radiation deduced from FUSE spectra (Harper et al. 2004). For further updated information 
on mass- loss mechanisms see, e.g., the review by Willson (2000). 

In the present work, we overcome these deficiencies by adopting a more physical picture. 
In our approach, the non-radiative energy input into the wind is assumed to be given by the 
turbulent energy density, within the chromosphere or underneath, possibly related to the 
manifestation of (magneto-)acoustic waves. This approach appears to be consistent with the 
major conclusion by Judge & Stencel (1991), who presented a detailed empirical analysis 
of the global thermodynamical properties of the outer atmospheres and winds of a set of 
well-studied cool giant and supergiant stars. They concluded that "[...] mass loss rates are 
not strongly dependent on the actual physical processes driving the winds [suggesting] that 
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nonlinear processes act to regulate wind energy fluxes". Furthermore, we assume that the 
mass loss rate depends on the characteristic chromospheric height, which dictates the amount 
of energy required to lift the wind out of the potential well of the star. This simplified model 
results in a mass loss formula akin to that by Reimers (1975, 1977). However, it contains 
two additional factors, one depending on the effective temperature and the second on the 
surface gravity of the star. This improves the agreement with observed mass loss rates for 
different types of stars, without the need to adjust the fitting parameter. In §2, we describe 
our theoretical approach. In §3, we discuss tests and applications of our new formula, and 
in §4, we give our conclusions. 



2. Theoretical Approach 

In our simplified model, we regard the wind to result from a spill-over of the extended, 
highly turbulent giant chromosphere and its reservoir of mechanical energy, possibly asso- 
ciated with waves. Even though no details of any theoretical models are considered, the 
following derivation will be guided by the assumption of Alfven waves, owing to their suc- 
cess in describing stellar winds as obtained for a Boo (K1.5 III) (Hartmann & MacGregor 
1980), C Aur A (K4 lb) (Kuin & Ahmad 1989), and a Ori (M2 lab) (Hartmann & Avrett 
1984; Airapetian et al. 2000). The relevant mechanical energy flux Fu may thus be due to 
magnetic energy generation, a consequence of non-isotropic chromospheric turbulence, or a 
combination of both. 

Turbulence is a well-known feature of stellar photospheres (e.g., Gray 1992, and refer- 
ences therein, among more recent literature), and of cool star chromospheres. E.g., Carpenter 
(1996) obtained empirical constraints on the chromospheric macroturbulence and flow ve- 
locities for various K and M (super-)giants based on C II] from HST-GHRS spectra, ranging 
from 24 km s" 1 (a Tau; K5 III) to 35 km s" 1 (a Ori; M2 lab), which are in principle sufficient 
to overcome the gravitational potential of the star. Nevertheless, the chromospheric turbu- 
lent energy density relevant for the generation of winds is not exactly known, largely because 
of the difficulty of distinguishing between isotropic and non-isotropic turbulence. Examples 
of wave-driven wind models have also been given. For instance, Airapetian et al. (2000) 
proposed a time-dependent, 2.5-D Alfven wave wind model, resulting in a time-averaged 
mass loss rate commensurate with the recent semi-empirical chromosphere and wind model 
by Harper, Brown, & Lim (2001) based on NRAO VLA radio data. 

Aside from the amount of utilized mechanical energy flux, the stellar mass loss rate 
is also expected to depend on the characteristic chromospheric radius -Rchr, which dictates 
the amount of wind energy dE Wind needed by a mass element dM = Mdt to overcome the 
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gravitational potential of the star. For the wind energy balance, we thus obtain 

GM^Mdt „ , - , 
^Wi„d ^ — ^ oc F M -4nRldt (1) 

with M as mass loss rate, and i?*, M* as stellar radius and mass, respectively, Fm as 
mechanical energy flux, and G as gravitational constant 1 . 

A large body of literature has been devoted to describing convective turbulence of stellar 
atmospheres and the generation of waves as function of the fundamental stellar parameters. 
Stein (1981) studied the generation of waves by turbulent motions in stellar atmospheres 
largely based on analytic means. He found that the acoustic wave energy flux is given as 
Fyi oc Tgg 1 (monopole term), T^ A (dipole term), and (quadrupole term) (see repre- 

sentation by Ulmschneider [1989]), noting that the monopole term is most closely and the 
quadrupole term is least closely related to mass loss generation. Models by Bohn (1984) 
deduce as temperature dependence for the monopole, dipole, and quadrupole terms, T c g 15 , 
T c g 78 , and T c g- 85 , respectively, and for the combination of those terms T^ 75 . Analytic models 
for magnetic wave generation reveal a temperature dependence of T e g 5 for the Alfven mode 
and T e 8 ff for the magnetic modes combined (Musielak & Rosner 1988). 

More recent work has resulted in vast improvements of these models, particularly with 
respect to the models for the stellar convection zones and the description of the turbulent 
frequency spectra — see, e.g., Musielak (2004) for a recent review on those results. Un- 
fortunately, the authors usually refrain from giving a fitting formula for the dependence of 
the wave energy fluxes on the stellar effective temperatures, supposedly because of the large 
number of free parameters, particularly in magnetic models. Nevertheless, based on the 
fact that the solar wind (e.g., Ong et al. 1997) and massive stellar winds (e.g., Rosner et 
al. 1995; Airapetian et al. 2000) are likely to be accelerated by the momentum deposition 
by Alfven waves, we consider a dependence as Tjg to be most representative. There is 
also another feature inherent to Alfven waves motivating this choice: they are known to be 
essentially non-dissipative in the lower and middle chromosphere — see discussion by Char- 
bonneau & MacGregor (1995) and Boynton & Torkelsson (1996), pointing to the prevalence 
of low-amplitude waves in those stars — implying that the same T e g-exponent holds in the 
region of wave generation as well as where the onset of mass loss occurs. Incidently, the 



1 Historically, a similar approach has been undertaken by Fusi-Pecci & Renzini (1975), who assumed that 
cool star winds are due to acoustic waves. Later on, this assumption has been invalidated by detailed model 
simulations, which showed that acoustic waves fail to support significant mass loss as they dissipate most 
of their energy immediately beyond the stellar photosphere, and are thus incapable of meeting the potential 
energy requirement of stellar winds (e.g., Cuntz 1990; Sutmann & Cuntz 1995, and references therein). 
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same exponent is also found as temperature dependence of Mg II h, k emission in stars (both 
dwarfs and giants) of minimal activity, a likely indicator of the overall chromospheric energy 
density (Buchholz, Ulmschneider, & Cuntz 1998, see their Fig. 15). 

If the mechanical energy flux Fm utilized for generating stellar mass loss is assumed as 
Fm oc T e g 5 , the surface-integrated mechanical energy flux Lm can now be expressed as 

L M = Fm-^RI oc F M -L*/T e 4 ff oc L* • T c 3 ff 5 . (2) 

Next we consider the characteristic chromospheric radius -Rchr- For cool giants and super- 
giants, no well-defined boundary between the chromosphere and the wind exists. Hence, we 
use the sonic point of the average velocity field as reference. For the well-studied K super- 
giant ( Aur (with log ~ 0.8), -Rchr is found to be close to 2R* (Baade et al. 1996), and for 
general giants and supergiants, (-Rchr ~ R*)/R* is assumed to vary as g' 1 , which gives 

R Chr = R*.(l + — ^ ) . (3) 

V 4300-0*/ V ; 

With the above temperature dependence of the mechanical energy flux (eq. [2]) and chro- 
mospheric radius i?chr (eq. [3]), we finally obtain as mass loss rate M, see eq. (1), 

M = V .^.(-^-f 5 + (4) 
' M* V4000 K/ V 4300 ■ gj V ; 

with R*, M*, and as stellar radius, mass, and luminosity given in solar units, and g* and 
g Q as stellar and solar surface gravity, respectively. This is, apart from the two new factors, 
indeed the old Reimers law. To satisfy the well-constrained RGB mass-loss of globular cluster 
stars (see § 3), the fitting parameter rj will be set to 8(±1) x 10 _14 M o yr _1 . 



3. Tests and Applications 

For the newly developed mass loss formula various tests and applications have been 
devised. In particular, we want to obtain insight into the importance of the new factors 
given by the stellar effective temperature T e s and gravity g*. In fact, for ordinary giants 
the two new factors do not make much difference, which explains the long-lasting success of 
the Reimers relation. In particular, the T e g 5 factor is, despite its high power, restricted in 
its impact by the small band of relevant effective temperatures (3000 to 4500 K), and the 
^-sensitive factor remains of the order of one for all but the smallest gravities. In fact, as 
previously discussed, the T e s exponent in eq. (4) is somewhat uncertain. However, due to 
the narrow band of relevant effective temperatures, the overall results would still stand if 
T„ 3 ff or were used instead. 
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Nevertheless, we have obtained evidence for the extra dependence on T eS from a compar- 
ative study of the RGB mass loss of globular cluster (GC) stars with very different metallicity 
and, accordingly, different effective temperatures on their RGBs. The mass lost on the RGB 
of a GC is very well constrained by modelling the stars at the blue end of the Horizontal 
Branch (HB), for which the HRD position is very mass sensitive. The remaining uncertainty 
is about 15 % in absolute terms and much better in relative terms. In fact, the long time (on 
a dynamic time-scale) spent on the RGB effectively evens out most of the inherent variability 
of these stellar winds. This is a big advantage over the kind of snapshots obtainable from 
observing individual winds directly. A residual star-to-star variation of < 20%, on the other 
hand, is sufficient to explain the full spread of an HB. 

As test cases, we consider two globular clusters, which are NGC 5904 and NGC 5927. 
NGC 5904 has a significant metal underabundance of [Fe/H] = -1.29 (normal for globular 
clusters, Z=0.001, see Fig. 1), whereas NGC 5927 only has a marginal underabundance of 
[Fe/H] = -0.37 (near Z=0.01, see Fig. 2). We use the photometric data and metallicities 
provided by Piotto et al. (2002) and have plotted our evolution tracks directly into their 
cluster CMDs. Clearly, the Reimers relation cannot reproduce both cases with the same t^r, 
while our new relation can! In particular, the extreme blue end of the Horizontal Branch 
(HB) of NGC 5904 (Fig. 1, bottom panel) demands HB stellar masses M HB ^ O.6OM , with 
a He-core mass of M c ~ O.49M . With an approximate initial mass of M, ~ O.86M , 
consistent with an age of about 12 billion years, this corresponds to a total RGB mass loss of 
0.26 M Q for the individual GC stars. With our new mass loss relation with rj = 0.8 x 10~ 13 , 
our evolution models achieve exactly this HB mass. To get the same result with the old 
Reimers law, we would need a rjn = 2.4 x 10~ 13 . On the other hand, the small extent of 
the HB of NGC 5927 (Fig. 2, bottom panel) demands HB stellar masses M HB ^ O.71M , 
possibly slightly more, but certainly not less, with a He-core mass of M c ~ O.48M . With 
Mi ~ O.99M , consistent with an age of about 11 billion yrs, this corresponds to a total RGB 
mass loss of 0.28 M — exactly as achieved by our evolution models with r] = 0.8 x 10~ 13 
and our new mass-loss relation. To get the same HB mass with the Reimers relation, we 
would need r/ R = 2.0 x 10~ 13 . On the other hand, a value of r/ R of 2.4 x 10~ 13 can clearly be 
ruled out: it would produce HB stars with M H b — O.63M , which are already far too blue. 

The role of the new, gravity-related factor appearing in our new mass loss relation can 
be explored by assessing luminous low-gravity stars, which play a very important role in the 
overall mass loss during the stellar lifetime. A good test candidate is the well-studied star 
a Ori, which has some circumstellar dust but is still below its critical luminosity for possessing 
a truly dust-driven wind. The mass loss rate was observed as M ^ s = 3.1 (±1.3) x 1O~ 6 M yr _1 
(Harper et al. 2001), assuming a distance of d = 131 pc (±25%) according to the Hipparcos 
parallax. The corresponding luminosity is = 5.4 x 1O 4 L , and with an angular diameter 
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of 56 mas, T eff is given as 3140 K, for which matching evolution tracks imply a mass of M* = 
1OM (±30% depending on L*). For these parameters, our new mass loss relation (with 
T) = 0.8 x 10~ 13 ) yields 2.2 x 1O~ 6 M yr _1 , while the Reimers law gives only 0.8 x 1O~ 6 M yr _1 
(with ?7r = 2 x 10~ 13 ). The distance uncertainty for a Ori affects all three mass loss rates 
(including M Q b s oc d 2 ), but not so much their ratios. For example, if d was 25% larger and 
M* ~ 13M , then M obs = 4.8 x 1O~ 6 M yr" 1 , our value would be 4.0 x 1O" 6 M yr" 1 , 
whereas the Reimers law would give 1.3 x 1O~ 6 M yr _1 . In any case, for such low gravities 
(logg ~ —0.35), Rchr/R* increases significantly, resulting in an extra boost of mass loss. 

4. Conclusions 

We derived a new semi-empirical relation for the mass loss of cool winds, which so far 
has frequently been described by Reimers' law. Physically, the Reimers relation appears to 
suggest a picture in which the wind material is lifted from the photosphere by using a certain 
fraction of the stellar luminosity — even though it is well known that, with the exception 
of molecule-driven and dust-driven winds, cool star winds are not related to any type of 
radiation pressure. This apparent contradiction has now been resolved. 

The new relation is based on theoretical arguments assuming that the wind results from 
the turbulent energy density, within the chromosphere or underneath, possibly related to 
the manifestation of magnetoacoustic waves as, e.g., Alfven waves. Furthermore, the mass 
loss rate is assumed to depend on the chromospheric extent, which dictates the amount of 
energy required to lift the wind out of the potential well of the star. A more detailed analysis 
shows that the new mass loss formula is not applicable to molecule- driven, dust-driven, and 
pulsational winds, as in those cases highly temperature-sensitive feed-back mechanisms exist, 
resulting in a steeper dependence of the mass loss rate on the stellar parameters, including 
the metallicity, which is not reflected by the new formula. Moreover, pulsational winds are 
more episodic in nature, whereas the new mass loss formula only descibes time-averaged mass 
loss behavior. Also note that the new mass-loss formula is not valid for stars like the Sun, 
where information exists that different types of mass loss processes, resulting in slow and 
fast wind, operate on different horizontal and vertical scales, which is beyond the theoretical 
framework of this paper. 

The new relationship mostly reproduces the original Reimers law, except that it includes 
two additional factors, which further improve the agreement with observed mass loss rates, 
especially for (super-)giants with very low gravity. This improved agreement can be inter- 
preted as an indirect validation of Alfven waves as being primarily responsible for mass loss 
generation in those stars. A highly sensitive calibration of the new relation's fitting factor 
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rj has been achieved by modelling the mass lost on the RGB by stars near the blue end of 
the Horizontal Branch, using two globular clusters of very different metallicity (NGC 5904, 
NGC 5927). Further studies, considering sets of well-studied K and M-type giant and super- 
giant stars, including comparisons with the various mass loss formulas from the literature, 
will be given in the near future. 

We are grateful to Z. E. Musielak, D. Schonberner, R. C. Smith, and B. E. J. Pagel for 
helpful discussions. The paper also benefited from comments by three anonymous referees. 
This work has been supported by NSF under grant ATM-0087184 (M. C). 
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Fig. 1. — Top panel: Color-Magnitude Diagram of the globular cluster NGC 5904. Bottom 
panel: Position of the Horizontal Branch with Z=0.001. 
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Fig. 2. — Top panel: Color-Magnitude Diagram of the globular cluster NGC 5927. Bottom 
panel: Position of the Horizontal Branch with Z=0.009. 



